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Abstract 

We consider domain walls generated through a cosmological phase transition, 
which interact non-gravitationally with light neutrinos (m v ~ 0(l)eV). At a redshift 
z > 10 4 the network grows rapidly and is virtually decoupled from the matter. As the 
friction with the matter becomes dominant, a comoving network scale close to that 
of the comoving horizon scale at z ~ lO 4 gets frozen in. During the later phases, the 
walls produce matter wakes of thickness d ~ 10/i -1 Mpc, that may become seeds for 
the formation of the large scale structure (30 — l30h~ l Mpc) observed in the Universe. 
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The discoveries of supergalactic structures at a scale comparable to lOO/i l Mpc 
and of an intergalactic medium already ionized at z ~ 5 prompted new interest in 
models of galaxy formation alternative to the standard cold dark matter (CDM) 
scenario. 1,2 As is well known, the CDM model finds it difficult to generate the very 
large clustering scales we observe, given the constraints from the microwave back- 
ground. This model also predicts that most of the galaxy formation should take 
place at a redshift 2 < 2, and can hardly explain the production of big quantities of 
ionizing radiation at early times. An alternative class of models is based on the fact 
that cosmological phase transitions could give rise to highly concentrated pockets of 
energy, that may become seeds for baryonic infall and lead to galaxy and structure 
formation earlier than what the standard model postulates. 

Cosmological phase transitions can result from the spontaneous breaking of a sym- 
metry associated with “Higgs-like” scalar fields. The symmetry breaking is usually 
considered to take place at a temperature close to the electroweak or the grand- 
unification (GUT) energy scales. The phase transitions can give rise to topological 
defects. Depending on the choice of the scalar fields one may, in particular, gener- 
ate 1-d or 2-d objects like “strings” or “domain walls”, formed of a “false vacuum' 1 
phase. As many studies show, these solitons can be associated with galactic and 
supergalactic scales and generate fluctuations in the distribution of matter. Typi- 
cally, the defects move at relativistic speeds and create matter fluctuations by their 
gravitational action. 3 

In the first model dealing with cosmological domain walls, Zel’dovich et al.(1974) 4 
considered a phase transition taking place at the GUT scale. The main results of 
the paper were discouraging. The walls produced were far too heavy, giving rise to 
unphysically large distortions of the microwave background radiation. The idea was 
therefore abandoned until a possible mechanism was found to form domain walls at 
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a much lower transition temperature. 5 It was found 6 that the surface energy density 
<7 of the walls is compatible with the constraints of the microwave background, by re- 
quiring a < 10 MeV 3 , but still gravitationally significant, if a ~ lOMeV 3 . Numerical 
simulations 7,8 showed in detail that the typical scale associated with the network is 
always close to that of the horizon, provided that the walls freely stretch under their 
surface tension. 

In a recent publication, 9 I discussed a variation of the previous domain wall models 
that could possibly generate both the large scale structure and give rise to a very 
early galaxy formation (at z ~ 10 — 30). A simple non-gravitational Lagrangian 
coupling was introduced between the scalar field associated with the walls and a 
component of the dark matter. We assumed that the domain walls gave rise to a 
symmetric potential barrier for the dark matter, resulting in elastic scattering. 10 As a 
simplifying assumption, the barrier was taken to be high enough to reflect elastically 
all the incident particles, regardless of their impact energy. The result of such an 
interaction is to keep the velocity of the walls to very non relativistic values (typically 
v ~ 10 -2 — 10 -3 c). Therefore, the network just stretches with the universal expansion, 
the comoving scale of the system being frozen. Summarizing the results, for this to 
occur, ^twaiis/^DM — 10~ 3 . The walls, by pushing ahead the dark matter, create a 
void region behind them and a high density wake in front of them, whose thickness 
can be of order d ~ lOh' 1 M pc. 

In the present paper, we modify the assumptions of the previous work. We consider 
the effects of a finite barrier height E a for the particle reflection. 10 We assume that 
particles hitting the domain wall with a relative kinetic energy E < E a are scattered, 
and for E > E a they pass through the barrier. As we show in ref. 10, this is always a 
good approximation, in the limit E 0 m, where m is the mass of the dark matter 
particles. We are going to limit our calculations to such a limit. In considering the 
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dark matter as a fermion gas in thermal equilibrium, we will find a link between the 
comoving scale of the network f and the mass m of the particles. 

We begin by deriving the equation of motion of the w T alls, under the assumption 
that they move through a homogeneous medium formed by a fermion gas in thermal 
equilibrium. The action of the dark matter gives rise to a pressure Pj on the domain 
walls. Deriving Pj is the first step of our calculation. 

If we define as the x-axis the normal to the moving surface, the momentum ex- 
change between wall and incident particles is A p ~ 2m^{v— v x ) (notice that, for the in- 
teraction to occur, (v -v x ) < 1, since E 0 < m), where v x is the normal component of 
the particle thermal speed in the background frame and 7 = (l-u 2 )~ 1/2 ~ (l-u 2 ) -1/2 . 
The hit rate is given by n(v — u r ), where n is the particle number density. It is con- 
venient to introduce the dimensionless variables y = m^v/T and y x = p x /T (where 
p x is the x- momentum of the particles in the rest frame and T is the temperature). 

We now introduce the thermal distribution f(y x ) of the particles in the momentum 
component p x (averaging over the other directions). 9 This function is defined so that 
I^oo f(Vx) d(y x ) = 1- Since the particles are fermions, 


,, x J_ f°° y± 

3C(3) Jo txp^Jyl 


Vx-dyL 


+ vi +1 


(1) 


where y L = p^/T (where p\ = p 2 - p\ and p is the particle momentum). Recalling 
that the reflections only occur if (v — v x ) -C 1, the pressure Pj is 


r rv /v+vo 

P/(y) = 2fT 2 n/m\ j (y - y x ) 2 f{y x )dy x - f (y - y x ) 2 f{y x )dy x , (2) 

Uy-y 0 J V 

where y 0 is a limiting value simply related to E 0 . In fact, in the rest frame of the 
wall, the maximum value of the incident particle momentum giving rise to reflection 
is approximately p 0 = yj2mE 0 (since E 0 < m). By boosting such value back to the 
universal rest frame we find y 0 = 7>/2 rnEl/T to a good approximation. 
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The first integral in eq.(2) is the momentum exchanged per unit time and area 
between the wall and the particles in front of it, while the second integral refers to 
particles hitting the wall from the back. Pj can be rewritten, by changing the variable 
of integration, as 


Pj = 27 T 2 n/m f y 2 f{y - y x )dy x - f y x f{y + y x )dy x 
0 J 0 


(3) 


by introducing y x = y — y x in the first and y x = y x —ym the second integral of 
eq.( 2 ). 

The range 7 y/2mE~ 0 T (y 0 > 1 ) corresponds to the high barrier limit, that we 

already calculated in ref.9. In this limit, Pj can be written as follows: 


12mm; 2 

mv T ,E 0 > mu 2 / 2 

(4) 

-vT* 

mv < T ,E 0 > mv 2 j2 

(5) 

T 


0 

mv ^>T ,E 0 < mv 2 / 2, 

(6) 


where v is the physical speed of the wall through the medium. 

We now want to study the limit 'y\/2mE 0 <C T (y 0 1), since we are interested 

in the behavior of the network at high z. It will turn out that the friction is ineffective 
at very early times, so that the network initially evolves only subject to its surface 
tension. Expanding / in power series around y, if 7 \j2rnEl/T <1 we get 


Pj = -47 nT 2 f'{y)/m F y\dy x . 

L J 0 


(7) 


The function f{y) is given by: 

_ 1 r 00 egpyV + 

3C(3) Jo Jyi + [earp^/y 2 + y 2 + l ] 2 

Substituting eq.( 8 ) into eq.(7) gives 


P f = - 47 nT 2 /mf'(y) f y*dy x = 

Jo 


2 m 7 E 2 v i 4 
7T 2 e y + 1 
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( 8 ) 

(9) 
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For y > 1 (recall that y = mv~f/T), the pressure gets exponentially small. Physically, 
when r > v tkermal » (2^,/m) 1 / 2 , the walls behave like virtually decoupled from the 

matter (see eq.9). 

In the limit y < 1, the T-independence of eq.(9) comes from the fact that the 
number of particles reflected by the walls is constant during the expansion, while the 

momentum exchanged is A p 0 ~ (Eo/m) 1 ^ 2 . 

The equation of motion of the domain walls, describing an infinitesimal wall seg- 
ment of curvature 1 /R (in physical coordinates), is given by 


« . ^ u- t 

7 2 u + 3-u H 

a (77 



( 10 ) 


where a is the surface energy density of the kinks. The derivatives are taken respect 
to the universal time t. This equation is valid if the wall thickness A <C R, which 
becomes valid soon after the phase transition. Eq.(10) is just the relativistic gener- 
alization of Newton’s second law divided by 7 < 7 , where a is the energy density of the 
walls. The damping term 3(a/a)t; derives from the universal expansion, while the 
r.h.s. is the tension that drives the motion of the walls. 

It is useful to recall the results of the limit Pf — ► 0, that has been already studied 
in several works. 4 ' 7 ' 8 It was shown that the network reaches a scaling regime of growth, 
the scale of the system being given by R = 2/3t a /3, with 0 = const. ~ 0(1), and 
Q ~ 1. A very small percentage of the network reaches highly relativistic speeds, and 
the average physical speed is constant. 

During the early stage of the wall evolution, there could be a phase during which 
T \J2mE 0 is satisfied. In this regime, Pf < nPE^v/ir 2 for most of the network 
(7 ~ 1). The cosmological damping term, instead, is Zav/a ~ t~ 1 v. Consequently, 
if the transition temperature T c \J2mE 0 and t~ l > m 2 El/x 2 (where t c is the age 
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of the Universe at the phase transition), the initial phase of the network evolution is 
frictionless and rapid. Under additional conditions that we will discuss, the comoving 
scale of the network freezes in at the value it had at the time t ~ (m 2 E 2 / 7r 2 ) j . 

Solving the general form of the equation of motion in the regime T ^/ < ltnE 0 
is a difficult task, since the expression for Pj is quite involved. There is, however, a 
wide range of parameters such that eq.(10) can be simplified. Given the assumption 
that 7 mv < T we can write the approximate expression: 


7 2 t> + 3 — v + 


m 


■Elv 


1 

R 


( 11 ) 


a Tr 2 a 

This equation can be studied analytically, provided that virtually all of the network 
obeys it. The linear dependence of Pj in v ensures a selfsimilar evolution throughout 
the process. 

With the onset of the friction, the walls slow down, until they stop evolving and 
just stretch conformally with the expansion. If Pj becomes dominant at z ~ 10 4 or so 
for all of the network, the freeze in gives a comoving scale of the order 10 — 10 2 Mpc, 
comparable to the horizon scale at that time. 

Since the network evolves self-similarly, we can perform an averaging procedure to 
transform eq.(ll) into an equation describing the evolution of the average ”interwalT 
distance. Such a procedure is analogous to the one used when dealing with domain 
walls in condensed matter physics. 11 

Multiplying eq.(ll) by v and averaging over the surface of the network we get 
/3a m 2 E 2 \ 2 1 2 . 1 ,v 

During the period in which Pj is negligible, (7 2 u • v) ~ (v • v) = d((v 2 ))/dt = 0. 
When Pj dominates, since v/v ~ —d/a, we get (v • v) ~ — a/a(v 2 ). This term can be 
always safely neglected. 
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It is convenient to express the parameters in comoving coordinates (R = af and 
v = ar ). From simple geometry, it follows 11 that (v/ R) = (r / r) = r 2 d(f 2 )/dt, where 
f is defined as the average comoving ’MnterwalF distance, related to the average 
curvature of the walls by the relation (djr = (l/r 2 ) 1/2 . Therefore we get, 

dr 1 \ a m 2 El 

—77 3 — I 

dt [ a a 

where (3 = 3 if there is one wall per horizon during the uncoupled phase. In what 
follows we will suppose that the dark matter is constituted by only one particle species. 
For an open Universe, the solution to eq.(12) can be written in the approximate form 

f = hh- x p^ln'l\\ + Ktsj/2)Mpc, (13) 

where we introduced K = ( m/\0eV) 2 {E o /10~ 4 eV) 2 ((T/lMeV 3 ) _1 for convenience, 
since typically K ~ 1 (see fig.(l)). The reason for the choice of the normalization of 
K respect to E a and a will be clear soon. The quantity t$f is the time of equivalence 
of matter and radiation: t 5J = ( t equiv /t(z = 10 5 )). After equivalence, the evolution of 
the network virtually stops. The values of the parameters involved indicate that the 
largest scales are obtained for small m. This is consistent with an open Universe, if 
the only forms of matter are neutrinos and baryons. 

The assumptions we have made in calculating f are selfconsistent if the network 
remains always confined to a range of speeds where Pf ~ v. This is approximately 
true if the network remains confined to the region where Pj is not exponentially 
suppressed. 12 Assuming that the distribution of the curvatures is roughly Gaussian, 
all but a negligible fraction (1%) of the wall sections have a comoving curvature 
r~ l < 2.5r -1 , and the condition above translates to a bound on K, namely K < 
1.5(m/10eV') 2 . For K > 1.5(m/10eV) 2 our analytical approach is invalid, since part 
of the network is decoupled and part is strongly coupled. In this complicated regime, 
a numerical simulation would be needed to determine the evolution of the network. 
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Considerations about the later evolution of the network can help us in further 
constraining the parameters. Up to this point we have calculated the evolution of 
the walls at T > ( mE 0 )^ 2 . As soon as the Universe cools down below T ~ (m£ 0 )U 2 , 
eq.(4) holds. Pj decreases rapidly, due to its steep T dependence, and v starts in- 
creasing. Shortly, v becomes bigger than the thermal velocity of the neutrinos and the 
walls are in the regime described by eq.(5). During this phase, as shown preliminarly 
in ref.9, the walls deplete the volume they sweep of dark matter, generating a matter 
wake in front of them. The speed of the walls increases at a slower rate: v ~ a. 

The details of the evolution through these different regimes are being investigated 
numerically. 13 Some order of magnitude estimate is already possible. The wails de- 
couple again as soon as v > (2E 0 /m) 1 ^ 2 . By that time (define it as tj), the comoving 
distance the walls have swept is roughly d/t u ~ v^td/dj = Vj, = (2 mE 0 ) 1 ^ 2 (where 
t u is the present age of the Universe, v & = v(tj) and aj = a(fj)). If we consider 
m = 1 — lOeU and d = 5 — 20 h~ l Mpc, we get E 0 = 10 -5-5 — 10 -3-5 eV\ From the con- 
straint on K , we also obtain a < (E o /10 _4 eV r )/1.5MeV r3 for our analytical approach 
to hold. The wall “re-decoupling” typically takes place at a<i ~ C?(10 -1 ). Since the 
network rapidly reaches relativistic speeds, we expect one or few wails to be within 
our present horizon. However, the wakes generated at t < tj should have gravitation- 
ally collapsed, in that involving also the baryon component. This collapse may been 
responsible for the formation of a first generation of protogalaxies at a z higher than 
the values predicted by the CDM model. 

The distortions of the cosmic microwave background radiation (CMBR) generated 
by walls and by the dark component of the wakes are beyond the present capability 
of detection. We give a fairly detailed analysis of the problem in ref.(9), and therefore 
we limit the present discussion to a brief summary of the main results. 

The biggest source of distortion associated wth walls and neutrinos is the grav- 
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itational potential of the matter wakes generated by the wall motion. Because of 
the Reese-Sciama effect, 14 these wakes give rise to a microwave background temper- 
ature distortion A T/T ~ Gp^R^X 1 ' 2 ~ 1°' 7 - An additional effect is due to the 
domain walls themselves, through the same mechanism. Since today there should 
be roughly one stretched and relativistic infinite wall within our horizon (plus col- 
lapsing wall bubbles), we expect A T/T ~ G<rt u ~ 10 -6 , for a ~ lA/eu 3 , as shown 
in previous work. 6 Both the effects are lower than the current limit A T/T ~ 10 -5 . 
The effects on the CMBR due to the baryons have not yet been treated. We expect 
that a major contribution to the distortions is due to the Zeldovich-Sunyaev effect, 
which is the scattering of low energy photons by free and hot electrons. This kind of 
CMBR distortion is typical of any model dealing with early formation of protogalactic 
objects. 

Special thanks to E. Kolb, A. Chakravorty, D. Schramm and A. Kosowsky, for the 
time devoted in reading the manuscript and their constructive remarks. I also thank 
N. Vittorio for remarks on the CMBR distortions. 
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Figure caption 


The comoving scale of the network f in h~ l M pc units, as a function of a. The 
friction begins to slow down the network at a ~ 10 -5 . The curves refer to different 
values of the parameter K and the neutrino mass m. For each value of m we take 
the limiting value K = 1.5(m/10eV) 2 . We also take 0 = 3 (see eq.(13)). 
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